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Forward calorimetry in the PHOBOS detector has been used to study charged hadron production 
in d+Au, p+Au and n+Au collisions at yf Snn = 200 GeV. The forward proton calorimeter detectors 
are described and a procedure for determining collision centrality with these detectors is detailed. 
The deposition of energy by deuteron spectator nucleons in the forward calorimeters is used to 
identify p+Au and n+Au collisions in the data. A weighted combination of the yield of p+Au 
and n+Au is constructed to build a reference for Au+Au collisions that better matches the isospin 
composition of the gold nucleus. The pt and centrality dependence of the yield of this improved 
reference system is found to match that of d+Au. The shape of the charged particle transverse 
momentum distribution is observed to extrapolate smoothly from p+p to central d+Au as a function 
of the charged particle pseudorapidity density. The asymmetry of positively- and negatively-charged 
hadron production in p+Au is compared to that of n+Au. No significant asymmetry is observed 
at mid-rapidity. These studies augment recent results from experiments at the LHC and RHIC 
facilities to give a more complete description of particle production in p+A and d+A collisions, 
essential for the understanding the medium produced in high energy nucleus-nucleus collisions. 

PACS numbers: 25.75.Dw, 25.75.Gz 


I. INTRODUCTION 

The PHOBOS detector [Tj at the Relativistic Heavy 
Ion Collider (RHIC) [2] is one of several experiments [3] 
[5] that have measured the invariant yield of charged 
hadrons in collisions of deuterons with gold nuclei at 
a nucleon-nucleon center of mass energy of = 

200 GeV. In the referenced papers, charged hadron pro¬ 
duction is studied as a function of both transverse mo¬ 
mentum (px) and collision centrality (a measure corre¬ 
lated with the impact parameter of the deuteron). The 
particle yields for px above about 1.5-2.0 GeV/c are sim¬ 
ilar to, or possibly slightly enhanced above, those ob¬ 
served in p+p collisions at the same energy [6], some¬ 
what reminiscent of the so-called Cronin effect seen in 


proton-nucleus collisions [7]. Previous analyses of the 
d+Au charged hadron spectra by PHOBOS [2] and the 
other RHIC experiments mm have demonstrated that 
this enhancement stands in stark contrast to the observed 
suppression of high px hadrons in the (central) Au+Au 
collision system at ,/s NN = 200 GeV [T3HI5] . Since no 
suppression is found in d+Au collisions, the effect seen 
in central Au+Au interactions has been interpreted as 
evidence of final state effects, in particular parton en¬ 
ergy loss. It should be noted that evidence of possible 
collective effects in systems such has d+Au and p+Pb 
have been found recently, but only for events with very 
high final state particle multiplicity (see, as one example, 
Ref. nn ). 

The choice of the reference system used in compar- 
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ing to Au+Au data, and of the centrality measure, are 
both of critical importance to the understanding of the 
observed suppression. The data and Monte Carlo (MC) 
simulations presented in this paper are used to study the 
choice of centrality measure, as well as the choice of ref¬ 
erence system. Centrality measures based on the multi¬ 
plicity of particles in the high-pseudorapidity region as 
well as on the number of spectators in the gold nucleus 
are examined. To study the chosen reference system, 
a calorimetry-based technique is used to identify, on an 
event-by-event basis, the subsets of d+Au collisions in 
which only the proton or only the neutron participated 
in the collision. Specifically, a calorimeter on the side of 
the interaction region where the Au beam exits is used as 
part of the determination of collision centrality while a 
second calorimeter on the other side is used in the selec¬ 
tion of n+Au and p+Au interactions. Similar tagging of 
the nucleon+Au component of the d+Au data has also 
been investigated by the PHENIX collaboration (T5J [lU]. 
These nucleon-nucleus collisions are used to construct 
an ideal reference system for comparison with Au+Au 
collisions. Further, the charged hadron yields of n+Au 
and p+Au are compared in order to study the ability 
of nucleon-nucleus collisions to transport charge to the 
mid-rapidity region. 



Figure 1. Top: Schematic overview of the PHOBOS Au- 
PCAL, also showing the shielding and ZDC. The solid 
(dashed) lines show the approximate trajectories followed by 
spectator protons from the Au nucleus with momenta of 100 
(50) GeV/c as they are bent by the DX magnet into the 
calorimeter. Bottom: Detailed expanded view of the PCAL 
region (right half of the upper figure), including the DX and 
DO accelerator magnets. The scale shown applies only to this 
detailed view. The shielding and ZDC detector are not shown 
in the bottom image. 


II. THE PHOBOS DETECTOR 


The PHOBOS experiment makes use of multiple detec¬ 
tor components to measure particles produced by colli¬ 
sions at RHIC. Silicon pad detectors near the interaction 
point are used for particle tracking and collision vertex 
determination, see Sect.[lV] Additional silicon pad detec¬ 
tors provide full azimuth and large pseudorapidity cov¬ 
erage, as described in Sect. |III| Collision triggering is 
provided by plastic scintillator arrays at high pseudora¬ 


pidity, see Sect. IIP and by calorimeters measuring the 
number of neutral spectator nucleons, described below. 
More detail on these subsystems may be found in Ref. [T]. 

To study nucleon-nucleus collisions, two calorimeters 
were added to the PHOBOS experiment prior to the 
2003 d+Au physics run at RHIC. These detectors extend 
the measurement of forward-going nuclear fragments. 
Complementing the pre-existing zero-degree calorime¬ 
ters (“ZDCs”) that collect energy from spectator neu¬ 
trons [20] , the proton calorimeter (“PCAL”) detectors 
measure energy from free spectator protons. Each PCAL 
detector is assembled from lead-scintillator bricks origi¬ 
nally constructed for the E864 experiment [21] at the 
AGS. The bricks are 117 cm in length with a 10 x 10 cm 
cross section facing the interaction point. Each brick has 
an array of 47 x 47 scintillating fibers running along its 


entire length. All of the fibers from a single detector ele¬ 
ment are read out by a Phillips XP 2262B phototube at 
the back. 

The PCAL detector on the Au-exit side of the collision 
(see plan view in Fig. I]) consists of an array 8 bricks wide 
by 10 bricks high. The d-exit side PCAL (not shown in 
the figure) is a small 2x2 array. As mentioned above, 
the former is used for centrality determination while the 
latter is used, along with the ZDC, for tagging n+Au and 
p+Au interactions. Both calorimeters are centered at the 
beam height and the smaller calorimeter is mounted with 
its elements at the same location transverse to the beam 
as the two closest elements shown in Fig. [l] 

Because of their higher charge to mass ratio (compared 
to the deuteron and Au nuclei, as well as nuclear frag¬ 
ments), spectator protons emerging from either side of 
the interaction are bent out of the beam pipe and into a 
PCAL detector by the RHIC DX-magnets. The primary 
purpose of these DX-magnets is to direct the deuteron 
and gold ion beams into and out of the interaction re¬ 
gion. 

The larger Au-PCAL covers a pseudorapidity region 
—3.6 < r] < —5.2 and therefore could be struck by pro¬ 
duced particles in addition to the spectator protons it was 
designed to detect. In order to prevent this, two shields 
consisting of 44 cm thick concrete blocks were installed 
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between the calorimeter and the interaction region. 

The energy coining from Au-side spectator protons 
(EPcal) is calculated using bricks in the the Au-PCAL 
which are located in the two rows at beam height, as 
well as the outer four columns away from the beam. The 
two rows at beam height are found to contain a major¬ 
ity of the hadronic shower energy in simulations of single 
nucleons having momenta comparable to nuclear thermal 
and fragmentation emission. The columns away from the 
beam supplement the shower containment. The remain¬ 
ing bricks, in columns near the beam but above and below 
beam height, are not included in EPcal. This reduces 
contamination from particles emitted in the neutron- 
induced hadronic showers which escape the ZDC. 

The Au-PCAL modules have been calibrated relative 
to each other using energy deposited by cosmic rays. Fast 
scintillator detectors are installed above and below the 
Au-PCAL detector to serve as cosmic ray triggers dur¬ 
ing dedicated calibration data taking. Modules in the 
d-PCAL have been calibrated relative to each other by 
minimizing the width of the single-proton peak in the 
d-PCAL energy distribution. 


III. COLLISION RECONSTRUCTION 

A. Collision Selection 


Deuteron-gold interactions are identified using a set 
of selection criteria designed to minimize background 
(i.e. beam-gas interactions) and enhance the sample of 
collisions which could produce particles inside the spec¬ 
trometer acceptance. First, at least one hit is required 
in both of the 16-scintillator arrays (see Ref. jTj for more 
details on this and other detector elements) which cover 
a pseudorapidity range of 3 < |?y| < 4.5. Then, the 
longitudinal collision vertex, as determined by a single¬ 
layer silicon detector covering the beam-pipe in the mid¬ 
rapidity range, is required to be within 10 cm of the nom¬ 
inal interaction point. Further, this vertex is required 
to be in reasonable agreement (within 25 cm) with that 
found by the simple timing difference of two sets of fast 
Cherenkov counters, located at —4.9 < ry < —4.4 and 
3.6 < 77 < 4.1 (ry > 0 being in the deuteron direction). 
Finally, events that appeared to have signals from ei¬ 
ther a previous or following collision are removed. If two 
events occur within 5 /xs, the later event is rejected as 
containing pile-up signals in the silicon. If two events 
occur within 500 ns, as determined using the fast trigger 
detectors, then both are rejected as pile-up. 



Figure 2. The event selection efficiency as a function of the 
ERing centrality variable. Grey points show the fraction of 
events simulated using AMPT that pass the event selection 
(see text for details). The black line represents a smooth fit 
to the points. 


B. Centrality Determination 

Two experimental observables have been used as cen¬ 
trality measures by the analysis presented in this paper. 
The first variable, “ERing”, is a measure of the total en¬ 
ergy recorded in “Rings”, endcap silicon detectors. The 
rings have nearly 27r coverage in azimuth and cover eta 
ranges of —5.4 < ?y < —3.0 and 3.0 < ry < 5.4. The 
second variable, EPcal, is described in Sect. [TlJ and mea¬ 
sures the energy of Au protons that do not participate 
in inelastic collisions with the deuteron. Thus, EPcal 
measures protons near beam rapidity, y = 5.36. 

The distribution of each of these variables in the d+Au 
data can be used to determine the fractional cross sec¬ 
tion centrality bins. Details on this procedure may be 
found in Refs. [551 [23]. The extraction of average values 
of collision parameters, such as the number of partici¬ 
pant nucleons (A part ), as well as the determination of the 
centrality-dependent efficiency of the collision event se¬ 
lection requires a set of simulations. Models of d+Au col¬ 
lisions from both the HIJING [55] and AMPT [55] pack¬ 
ages have been studied. The detector simulation has been 
performed using the GEANT package [26]. In addition 
to lVp art , other centrality parameters have been studied 
using these simulations, including A^“ t and A part , the 
number of participants in the gold and deuteron, respec¬ 
tively, A co n, the number of nucleon-nucleon collisions in 
the interaction, and iq the average number of collisions 
per deuteron participant. 

The efficiency of the collision selection can be deter¬ 
mined from the simulations as a function of the chosen 
centrality variable. This is done by counting the frac¬ 
tion of simulated events that pass the event selection as 
a function of centrality. Because the event selection con¬ 
tains a vertex cut, the fraction is calculated as the num- 
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ber of events passing the event selection divided by the 
number of events having a true interaction vertex within 
10 cm of the nominal interaction point. The efficiency 
as a function of the ERing centrality variable, obtained 
using AMPT simulations, is shown in Fig. [2j Note that 
the efficiency does not approach unity, even for central 
events, due to the small acceptance of the detectors that 
determine the collision vertex as part of the trigger. 

This efficiency is used to unbias the centrality variable 
distribution measured in the data. The unbiased distri¬ 
bution is then divided into fractional cross section bins, 
using the method described in Ref. EH 

The efficiency function presented in Fig. [3] is also used 
to correct, on an event-by-event basis, the measurements 
of the charged hadron spectra presented in this article. 
This accounts for the variation of the selection efficiency 
within a centrality bin, whereas the application of the 
average efficiency in a centrality bin would not. 

For both the HIJING and AMPT collision generators, 
a Glauber model has been used to determine the average 
values of centrality parameters, such as -/Vp art , which can¬ 
not be measured directly. A Hulthen wave function [38] 
has been used to model the deuteron profile, while the 
gold nucleus density has been modeled using a Woods- 
Saxon distribution. The value of the inelastic nucleon- 
nucleon cross section used in the Glauber model is 41 mb. 
The average value of the chosen centrality parameter can 
then be determined for each fractional cross section bin; 
for details on this procedure, see Ref. [29] . 

The systematic uncertainties of the various (unbiased) 
centrality parameters, such as -/V par t, have been studied. 
The dependence on simulations has been quantified by 
varying the centrality efficiency, for example, that shown 
in Fig.[2]for ERing centrality bins. The amount by which 
the efficiency can vary is estimated by dividing the sim¬ 
ulated events into vertex bins. The dependence on the 
deuteron wave function has been studied by using both 
a Hulthen wave function, as well as a Woods-Saxon dis¬ 
tribution. The uncertainty of the centrality parameters 
resulting from the choice of collision simulation model has 
been studied by comparing to simple Glauber MCs. Un¬ 
certainties in using the efficiency function to unbias the 
centrality parameters have been accounted for by smear¬ 
ing the centrality measure (i.e. ERing) prior to applying 
the efficiency correction. Finally, the centrality parame¬ 
ters coming from different collision simulation packages 
are compared. 

The centrality parameters determined from ERing cen¬ 
trality bins are presented in Table [TJ The values for 
p+Au and n+Au tagged events, described in Sect. |IIID| 
are also shown. The systematic uncertainties of (-/Vp art ) 
and (N con ) are typically slightly different, with that for 
(iVcoii) usually larger. The table lists the larger of the 
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Figure 3. The correlation between EPcal and ERing used to 
obtain EPcal centrality bins. 

two uncertainties. 


C. Proton Calorimeter Centrality Determination 

The Au-PCAL detector facilitates the determination of 
the centrality of d+Au collisions using a variable, EPcal, 
which is independent of the measured multiplicity. As 
has been shown previously |221123] , multiplicity measure¬ 
ments in a particular region of pseudorapidity may be 
biased if the centrality of collisions is determined using 
(multiplicity based) observables in a similar pseudora¬ 
pidity region. The ERing observable is measured at high 
pseudorapidity, allowing measurements at nrid-rapidity 
to be minimally biased by such auto-correlations. Cen¬ 
trality derived from the number of spectator nucleons 
should be free of such biases. A measurement of the 
charged hadron spectral shape in centrality bins from 
both ERing and EPcal is presented in Sect. |VI| 

Centrality bins could, in principle, be derived from 
EPcal signals using the same procedure as for the other 
observables. However, the breakup of the gold nucleus 
is not modeled by either the MC event generators, HI¬ 
JING and AMPT, or the GEANT detector simulation. 
As a result, an alternative procedure has been developed 
that exploits the monotonic correlation in the d+Au data 
between the EPcal signal, and the signal of another (well- 
modeled) detector, ERing. This correlation is shown in 

Fig-i 

The method for deriving event selection efficiency for 
a given value of EPcal uses the known efficiency of ER¬ 
ing. Using each event in the data, two distributions of 
EPcal are generated: one simply counting events and one 
counting events but weighted by the inverse of the known 
efficiency of the correlated observable, 1 1 (-ERing- The ef- 
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Figure 4. The event selection efficiency as a function of the 
EPcal centrality variable. Points represent the ratio between 
the number of events in an EPcal bin and the number of events 
expected for a perfectly efficient detector, obtained using the 
ERing efficiency function (see text for details). The black line 
is a smooth fit to the points. 
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Figure 6 . fVpart dependence on ERing in the MC. The white 
line shows the fit to the mean N palt in each ERing bin. 



Figure 5. The centrality bins obtained using the known 

ERing efficiency. Each slice of the histogram shows the dis- Figure 7. The .Apart distribution in the 0-20% central EPcal 

tribution of EPcal within the specified fractional cross section bin found using the .Apart vs ERing fit method (open his- 

bin. The kinks in the shaded histograms arise from the edge togram) compared to that from the weighting method (grey 

of a fractional cross section bin falling inside a histogram bin. histogram). Each distribution is (independently) normalized. 


ficiency as a function of EPcal is determined from the 
ratio of the simple-count distribution divided by that us¬ 
ing weighted counts. This efficiency is used in the stan¬ 
dard procedure to evaluate EPcal cutoff values for the 
centrality bins. Figures [4] and [5] show the event selection 
efficiency as a function of EPcal and the resulting EPcal 
centrality bins, respectively, obtained by using ERing. 

Two different procedures have been developed to es¬ 
timate the average number of nucleons participating in 
the inelastic collision, 7V part , for a given EPcal central¬ 
ity bin. Both procedures exploit the correlation of EPcal 
with ERing and then of ERing with 7V part . The same pro¬ 
cedures are used to estimate other collision parameters 
as well, such as the number of nucleon-nucleon collisions, 
iVcoii, or the impact parameter, b. 

The simpler approach involves fitting the mean 7V part 


in small bins of ERing, as shown in Fig. [6] The fit is used 
to estimate the average value of fV part given the value of 
ERing in an event. These values are then used to obtain 
A r par t distributions for each EPcal centrality bin. 

The second approach begins by dividing the ERing dis¬ 
tribution for events in a given EPcal centrality bin by 
the distribution for all events in order to determine the 
probability of any particular value of ERing in that bin. 
Then, for each centrality bin, all MC events are weighted 
according to the appropriate probability for their value 
of ERing and the distribution of N pajIt is determined with 
these weights applied. 

The results of the two techniques are compared for the 
most central EPcal bin in Fig. [TJ In the first procedure 
(open histogram), the spread of A^ part in the resulting 
distribution depends only on the width of the correla- 
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tion of EPcal and ERing, while in the second (grey his¬ 
togram) it is also affected by the correlation of ERing and 
7Vp a rt ■ The latter is almost certainly an overestimate of 
the width of N part in a given centrality bin, while the 
former may underestimate the spread. However, in the 
analysis of spectra and yields, this difference in width 
is only significant to the degree that it affects the mean 
value. The differences of the means found using the two 
techniques are included in the systematic uncertainty es¬ 
timate for the values of N part . Analogous systematic un¬ 
certainties are determined for the other centrality param¬ 
eters, such as iV col i or b. The weighting and fit procedures 
differ by about 5% in central d+Au and about 25% in pe¬ 
ripheral p+Au. 

The systematic uncertainty inherent in the procedure 
used to determine centrality from the EPcal variable has 
been studied. This has been done by applying the indi¬ 
rect procedure described above for EPcal to well modeled 
detectors at mid-rapidity, for which the direct procedure 
described in Sect. Ill B| can also be used. Discrepancies 
between centrality parameters obtained via the direct 
and indirect methods are used to quantify the system¬ 
atic uncertainties on this procedure. This uncertainty is 
in addition to those described in Sect. mm 

The centrality parameters found in the EPcal central¬ 
ity bins are presented in Table [T] The parameters have 
been determined using the weighting method. The table 
also lists the values for p+Au and n+Au tagged events, 
which are described in Sect. mm 


discussed below. 

The subset of d+Au collisions matching p+Au and 
n+Au interactions have been identified through the ob¬ 
servation of the spectator nucleon of the deuteron. The 
deuteron spectators are measured in PHOBOS using 
both the PCAL and ZDC detectors on the deuteron-exit 
side of the collision. Qualitatively, a collision in which 
the d-PCAL registered a spectator and the d-ZDC did 
not is labeled an n+Au interaction (and vice-versa for 
p+Au interactions). 



Figure 8. Regions used to study the characteristics of events 
with different total charge deposited in the d-PCAL (color 
online). Region 0 is the black colored bin located at the lowest 
detected d-PCAL signal. 


D. Deuteron-Nucleon Tagging 

The low binding energy of the deuteron nucleus 
(1.11 MeV per nucleon) facilitates the analysis presented 
in this paper. Because the deuteron is so weakly bound, 
it is possible for the nucleons to be relatively far apart at 
the moment the deuteron collides with the gold nucleus. 
This can result in only one nucleon of the deuteron partic¬ 
ipating in the (inelastic hadronic) collision. Furthermore, 
the binding energy is so small compared to the beam en¬ 
ergy that the remaining spectator nucleon can emerge 
from the collision almost completely unperturbed. Thus, 
such a collision can be treated as an effective collision 
between a single nucleon and a gold nucleus. 

Although the size of a deuteron is relatively large, the 
proton-neutron separation is typically not larger than 
the size of the Au target. As a result, the nucleon- 
gold collisions that form a subset of the deuteron-gold 
data are not equivalent to minimum bias nucleon-gold 
data. Rather, they are biased towards more peripheral 
interactions. Further investigations of this bias and the 
techniques used to address it in the present analysis are 



ERing (a.u.) 


Figure 9. The ratio of the ERing distribution for events in 
each region of Fig. [8] to that of region 7, which is partially 
under the proton peak. The color of the lines follows the 
same scheme as used in Fig. [8] Region 0, in which no energy is 
deposited into the d-PCAL, shows a bias toward higher values 
of ERing, which is associated with more central collisions. 

The observation of a spectator by one of the detectors 
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Parameter 

System(s) 

0-20% 

ERing 

20-40% 

Bins 

40-70% 

70-100% 

EP( 

0-20% 

:al (from 
20-40% 

3Ring) E 
40-70% 

ins 

70-100% 

(b) (fm) 

d+Au 
p+Au, n+Au 

3.3(1.4) 

6.1(1.4) 

4.7(1.5) 

6.4(1.3) 

6.3(1.4) 
7.2(1.3) 

7.6(1.3) 
8.0(1.3) 

4.1(1.8) 

6.9(1.4) 

4.9(2.0) 
7.2(1.4) 

6.0(1.9) 
7.6(1.4) 

7.3(1.6) 

7.9(1.3) 

(-^Vpart) 

d+Au 
p+Au, n+Au 

15.4(3.8) 
9.4(3.4) 

10.6(2.9) 
7.7(2.5) 

6.3(2.4) 
4.7(1.9) 

3.1(1.3) 
2.7(1.0) 

12.8(4.9) 
5.8(3.0) 

10.4(4.9) 
4.8(2.6) 

7.4(4.3) 
3.9(2.1) 

4.1(2.5) 

3.0(1.4) 

(Aeon) 

d+Au 
p+Au, n+Au 

14.5(4.2) 
8.4(3.4) 

9.4(3.3) 

6.7(2.5) 

5.0(2.5) 
3.7(1.9) 

2.0(1.2) 

1.7(1.0) 

11.8(5.2) 
4.8(3.0) 

9.3(5.2) 

3.8(2.6) 

6.1(4.5) 
2.9(2.1) 

3.0(2.5) 

2.0(1.4) 

M 

d+Au 
p+Au, n+Au 

7.6(2.1) 

8.4(3.4) 

5.2(1.8) 

6.7(2.5) 

3.3(1.5) 
3.7(1.9) 

+7(0.9) 

1.7(1.0) 

6.3(2.6) 

4.8(3.0) 

5.2(2.5) 

3.8(2.6) 

3.7(2.3) 
2.9(2.1) 

2.2(1.4) 

2.0(1.4) 

Sys. Error 

d+Au 
p+Au, n+Au 

7.5% 

10% 

10% 

12% 

15% 

17% 

30% 

31% 

15% 

31% 

15% 

31% 

20% 

31% 

30% 

31% 


Table I. Centrality parameters determined using ERing- and EPcal-based centrality bins and AMPT collision simulations. 
Centrality bins represent the fraction of the total d+Au cross section, even for the p+Au and n+Au collision systems (see 
Sect. Ill El. Values in par enthe ses are the RMS of their respective parameters. For the EPcal bins, the weighted ERing method 
has been used (see Sect. Ill Cl. ( b) is the average impact parameter, (A part ) is the average number of participant nucleons, 
(fVcdi) is the average number of collisions and {v) is the average number of collisions per deuteron participant. The last row 
lists systematic uncertainties in (A r par t) and (A co n). See text for discussion. 


is established from the amount of energy deposited in 
that calorimeter. Because the response of these calorime¬ 
ters has not been simulated in the PHOBOS d+Au MC, 
the efficiency and purity of the chosen signal cuts can¬ 
not be studied directly. Instead, the effect of the cuts on 
an independent centrality measure (ERing) has been ex¬ 
plored. This alternative method is motivated by the ex¬ 
pectation that tagging nucleon-nucleus collisions should 
produce a data set that is biased toward interactions with 
larger impact parameters than the full d+Au data set. 

The distribution of energy deposited in the d-PCAL 
is shown in Fig. [8] which has been divided into an ar¬ 
bitrary set of regions numbered 0-8. While regions 7 
and 8 show evidence of a proton peak in the d-PCAL, 
events from all regions with non-zero energy deposition 
(regions 1-8) show similar centrality characteristics, as 
will be discussed below. Only events in region 0 show 
a bias toward more central collisions and are therefore 
assumed to completely lack a proton spectator. 

The presence or absence of a centrality bias in the re¬ 
gions displayed in Fig. [H] is seen in Fig. [9j which shows 
the variation in the shape of the ERing distribution for 
events depositing different amounts of charge in the d- 
PCAL. Each line represents the ratio between a particu¬ 
lar region of Fig. [8] and region 7. Collisions that deposit 
no energy into the d-PCAL show a striking bias towards 
more central (higher ERing) interactions. Collisions in 
regions 1-8 all show similar ERing distributions. This 
suggests that any amount of energy deposited into the 
d-PCAL is due to a proton spectator from the deuteron. 

Furthermore, the observation that the shape of the 
ERing distribution is the same for all collisions which 


deposit energy in the d-PCAL supports the idea that 
these collisions are all of the same type, namely n+Au. 
As expected, the centrality of d+Au and tagged n+Au 
collisions differ, but the centrality of n+Au does not de¬ 
pend on the amount of energy that the spectator proton 
deposits in the calorimeter. 



Figure 10. The d-PCAL signal versus the d-ZDC signal. The 
boxes (at high d-PCAL, low d-ZDC and vice-versa) show the 
regions in which collisions can be identified as p+ Au or n+Au. 
Note that the quadrant near (0,0) is not used to identify pure 
d+Au collisions, since it also contains nucleon-nucleus colli¬ 
sions (as the calorimeters are not perfectly efficient). 

A similar procedure has been followed in order to de¬ 
termine the range of energy deposited in the d-ZDC that 
corresponds to a neutron spectator. The final regions 
in which n+Au and p+Au interactions are identified is 


























shown in Fig. |10| The minimum value of energy deposi¬ 
tion in the d-PCAL is well above region 0, but ensures a 
very clean n+Au sample. 


E. Centrality of Nucleon-Nucleus Collisions 

The centrality of the tagged p+Au and n+Au collision 
data sets are quantified (by parameters such as lV co n) 
within the fractional cross section bins determined for 
d+Au. This is necessary because the forward calorime¬ 
ters are not included in the simulations of the detector 
response, which precludes an event tagging procedure 
based on the simulated energy deposition of those de¬ 
tectors. 

Within a d+Au fractional cross section bin in the MC, 
the centrality parameters of tagged events, such as -/V part 
in p+Au collisions, are obtained using the true subset 
of simulated d+Au events identified as p+Au or n+Au. 
These subsets are identified by the presence of a neutron 
or proton, respectively, emerging from the collision with 
a longitudinal momentum of 100 GeV/c. 

The use of a tagging procedure based on true simulated 
momenta is valid under the assumption that the event 
tagging procedure used in data has an efficiency that does 
not depend on the centrality of the collision. That is, 
the average value of N co \\ in p+Au is the same whether 
the tagging efficiency is 80% or 100 %, as long as the 
tagging procedure does not alter the shape of the N co \\ 
distribution (but merely scales its normalization). 

The validity of this assumption rests on three reason¬ 
able conjectures. First, that it is not possible for a nu¬ 
cleon to both interact inelastically and to still deposit a 
measurable amount of energy into a forward calorimeter. 
Note that the Au-PCAL acceptance covers only protons 
having no transverse momentum and a longitudinal mo¬ 
mentum |p z | > 20 GeV/c, and that the smaller d-PCAL 
will observe only protons with even higher momenta. Sec¬ 
ond, that if a spectator nucleon is present, it will be de¬ 
tected by a forward calorimeter with an efficiency that is 
independent of the collision centrality. Finally, that the 
forward calorimeter on the deuteron side detects only 
deuteron spectators and not produced particles. 

These conjectures imply that the tagging efficiency is 
independent of the centrality of the d+Au collisions. This 
allows nucleon-nucleus collisions to be extracted from 
the d+Au (AMPT) simulations, analogous to the tag¬ 
ging procedure used for data. Simulated nucleons emerg¬ 
ing from the interaction at the nucleon beam energy are 
taken to be spectators. The centrality parameters ex¬ 
tracted using this method are presented in Table jTj 

The p+Au and n+Au events from the simulations have 
been used to obtain a rough estimate of the tagging ef¬ 


ficiencies. Taking the ratio of p+Au (n+Au) to d+Au 
collisions that passed the event selection in the simula¬ 
tion gives the fraction of p+Au (n+Au) events in the 
d+Au sample that would be tagged with a perfectly effi¬ 
cient detector. Dividing the actual ratio of tagged p+Au 
(n+Au) to d+Au events found in the data by the frac¬ 
tion expected from simulation gives an estimate of the 
efficiency. It is found that ~ 63% of p+Au interactions 
and ~ 46% of n+Au interactions are tagged using the 
procedure described above. The lower n+Au efficiency 
may be at least partly due to the relatively large mini¬ 
mum d-P CAL e nergy required in the tagging procedure 
(see Sect. HID). 


IV. HADRON SPECTRA EXTRACTION 


The transverse momentum spectra of charged hadrons 
have been extracted from tracks reconstructed using hits 
in the 16 layers of silicon detectors that make up the two- 
arm magnetic spectrometer. Hit position information is 
obtained both inside and outside the 2 T magnetic field. 
Details of the vertex determination and particle track¬ 
ing, as used in previous PHOBOS d+Au hadron spec¬ 
tra analyses, have been described in Refs. mm- How¬ 
ever, the current studies make use of an expanded set 
of data and an updated reconstruction procedure. As 


the d+Au collision trigger (described in Sect. Ill A) does 
not include a high px particle trigger, as employed in [B], 
a less biased data sample has been used in the present 
analysis. To improve the efficiency of the particle recon¬ 
struction, the final minimization step of the tracking has 
been performed multiple times for each track. This helps 
prevent the reconstruction from falling into a local min¬ 
imum, which reduces the number of both poor-quality 
track fits as well as ghost tracks. 

In an effort to more accurately apply acceptance and 
efficiency corrections, several changes have been made 
to the procedure used to extract the hadron momen¬ 
tum spectra described in Ref. [B]. First, the geomet¬ 
rical acceptance and tracking efficiency correction have 
been applied separately for each of the two spectrometer 
arms. To account for acceptance effects as accurately as 
possible, the correction factors as a function of px have 
been applied as interpolated spline functions of the track¬ 
embedding results (described in Ref. [IB]), rather than 
as smooth analytic functions. Further, the minimum px 
of acceptable tracks has been lowered to correspond to 
the px value at which the acceptance and efficiency cor¬ 
rections are roughly 30% of their maximal value. This 
leads to a minimum px value of 0.3 to 0.4 GeV/c, de¬ 
pending on the longitudinal collision vertex position, for 
hadrons bending towards higher -77 (out of the PHOBOS 
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spectrometer acceptance) and a minimum px of about 
0.1 GeV/c for hadrons bending towards negative 77 . Cor¬ 
rections for dead and hot channels in the spectrometer 
have also been applied independently for each spectrome¬ 
ter arm, to account for discrepancies on the level of 1 % in 
the hot and dead channel fraction of the two arms. The 
number of ghost and secondary tracks passing the recon¬ 
struction cuts are corrected for as a function of px- Due 
to improvements in the reconstruction software since the 
publication of Ref. jb] , these corrections are on the order 
of 1%. Finally, corrections have been applied for the mo¬ 
mentum resolution of the tracking and the variable bin 
sizes of the spectra. These corrections are determined 
using a dedicated simulation of single particles through 
each spectrometer arm to determine the distribution of 
reconstructed transverse momentum in each (true) px 
bin. 


The efficiency of the event selection described in 
Sect. Ill A depends on centrality, particularly for periph¬ 
eral events. Spectra uncorrected for this effect would cor¬ 
respond to an ensemble of events with a biased (higher) 
number of participants, rather than to a minimum bias 
selection using the same centrality binning. Instead, 
the efficiency determined as a function of centrality (see 
Fig- § is used to correct the spectra. 


The spectra of charged hadrons for d+Au, n+Au and 
p+Au collisions are presented in Fig. El in four bins of 
d+Au centrality, as determined by the ERing variable. 
For n+Au and p+Au, the same ERing cuts were used as 
for d+Au. Therefore, these do not correspond to bins of 
the listed fractional cross-section for nucleon-gold inter¬ 
actions. Note that the difference in the px range between 
d+Au and the nucleon-nucleus spectra is simply due to 
fewer p+Au and n+Au collisions being collected com¬ 
pared to d+Au. 

Systematic uncertainties on the measured charged 
hadron spectra have been quantified using the data. The 
largest correction, the acceptance and efficiency of the 
tracking, is the source of the largest systematic error 
(about 8 % at px = 2 GeV/c). This error has been esti¬ 
mated by comparing the yield in different subsets of the 
data for which the particle spectrum is expected to be 
the same. For example, the charged hadron yield of data 
taken with the spectrometer magnet in the positive po¬ 
larity is compared to that of data taken with the magnet 
in the negative polarity. Similarly, yields measured sep¬ 
arately in each spectrometer arm have been compared in 
order to derive uncertainties arising from the dead and 
hot channel correction. With these corrections applied 
separately to each arm, the systematic uncertainty on 
this effect is reduced to < 3% from ~ 10% in the previ¬ 
ous analysis [5]. 


For corrections in which such studies are not possible, 


the uncertainties are taken to be of the same order as 
the corrections themselves. At px = 2 GeV/c, this gives 
a ghost track uncertainty of 1 %, an uncertainty on the 
effect of secondary tracks of 3% and an uncertainty on the 
momentum resolution and momentum binning correction 
(which are applied together) of about 3.5%. 

Uncertainty on the yield of nucleon-nucleus collisions 
due to tagging has been estimated. This is done by vary¬ 
ing the d-PCAL and d-ZDC cuts used to tag events, 
which is expected to impact the number of interactions 
in the data set, but not the yield of those interactions. 
The total systematic uncertainties for the charged hadron 
spectra are shown in Fig. [12] 

The charged hadron spectra are used to derive the mul¬ 
tiplicity near mid-rapidity for d+Au, p+Au and n+Au. 
Spectra are modeled with the following functional form 


1 d 2 N ch 
27rpx dpxd ?7 



+ 


B exp 



(1) 


In the actual fit, parameter A in Eq. [l]is replaced by its 
value in terms of the analytically integrated yield dN c h/d>] 
and the other four parameters. 


A = 


(n — l)(n — 2 ) 

2t rp§ 


dNch/dq- 

27ri?T(m w + T)e 


-m,/T 


( 2 ) 


This allows both the value of dN c h/dq and its statistical 
uncertainty to be obtained directly from the fit. Sys¬ 
tematic uncertainties on the multiplicity are obtained 
by simultaneously shifting each point in the spectra to 
the limit of its individual systematic error and extract¬ 
ing dN c h/dq. The resulting systematic uncertainty on the 
integrated yield is about 9%. 

The charged particle multiplicity near mid-rapidity, at 
( 77 ) = 0.8, is shown in Fig. 13 for d+Au, p+Au and n+Au 


as a function of -A part . The number of participants is de¬ 
termined using ERing centrality bins, since the ERing 
measurement of particles far from mid-rapidity has been 
shown to introduce at most a minimal bias on the mea¬ 
surement [30]. A consistent dependence of the multiplic¬ 
ity on ./Vp art is observed across all three collision systems. 


V. AN IMPROVED REFERENCE SYSTEM 

The yield of hadrons in d+Au collisions has played a 
vital role in the investigation of particle production in 
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d+Au n+Au p+Au 



Figure 11. The invariant yield of (h + + h~)/2, h + , and h“ in four centrality bins determined for d+Au using the ERing 
centrality variable. The spectra for d+Au, n+Au, and p+Au are shown in separate columns. Due to the use of identical 
ERing cuts in all cases, the different data sets do not correspond to the listed fraction of the total inelastic cross section for 
nucleon-gold interactions. See text for details. Only statistical errors are shown. The spectra are obtained using particles that 
have a pseudorapidity 0.2 < p < 1.4. 


high energy Au+Au collisions. The nuclear modification 
factor, i?Y, of a collision system, X, given by 

d . _ d 2 N x / dpT drj 

A {Ncow) d 2 N pP /dpTdr) 

where X =Au+Au, d+Au, etc., has been used to test the 
scaling of the high-px hadron yield with the number of 
binary nucleon interactions occuring during the collision. 
The nuclear modification factor of nucleus-nucleus colli¬ 
sions at RHIC has been studied extensively for Au+Au 


interactions at = 39 GeV [32] 62.4 GeV [321134] . 

130 GeV mm and 200 GeV QZU E2 ES[38], as well 
as for Cu+Cu interactions at ,/s NN = 22.4 GeV [35], 
62.4 GeV [35] and 200 GeV [35] Ej. 

One of the fundamental conclusions drawn from ex¬ 
amination of the nuclear modification factor is that the 
production of high-px charged hadrons in central Au+Au 
collisions at ys NN = 200 GeV is highly suppressed with 
respect to binary collision scaling [16] . However, it 
cannot be known from the nucleus-nucleus data alone 
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Figure 12. Contributions to the relative systematic uncer¬ 
tainty associated with uncertainty in the hadron spectra cor¬ 
rections. 



(N > 

part 


Figure 13. The measured dN c h/dr] at (rf) = 0.8 in different 
collision systems obtained using ERing centrality bins. Sys¬ 
tematic errors are shown as error bars for the d+Au measure¬ 
ments; statistical errors are negligible. Systematic errors on 
the nucleon-nucleus measurements are not shown, but are of 
similar order. Asterisk symbols show PHOBOS multiplicities 
at mid-rapidity from Ref. 13011311 . 


whether the suppression is due to initial [42] or final [43] 
state effects. Nucleon-nucleus collisions at the same cen¬ 
ter of mass energy would provide a control experiment 
capable of distinguishing between the two possibilities, 
as such collisions should provide a nucleus in the same 
initial state but should not produce an extended medium 
in the final state. At RHIC these studies have been 
performed using d+Au rather than nucleon-nucleus col¬ 
lisions (snamaiiara- The assumption was made that, 
due to the small size and weak binding of the deuteron 
nucleus, d+Au collisions would provide as good a control 


experiment for Au+Au interactions as nucleon-nucleus 
collisions. 

This assumption can be tested using tagged p+Au and 
n+Au collisions to construct an improved reference for 
Au+Au collisions. Previous studies performed by the 
NA49 collaboration [?5l @B] have suggested that hadron 
production of nucleus-nucleus collisions may be better 
understood through careful consideration of the neutron 
content of the nucleus. Taking into account the fact that 
a gold nucleus consists of 60% neutrons and 40% protons, 
an improved nuclear modification factor for comparison 
to Au+Au can be defined as: 


R 


pnAu 


ni W a 7^)/(^qii) m " 

dN^/dn 

(dA^V^/jiVcoii) 

dNg/dn 


nAu 


(4) 


where (N co n) pAu is the average number of collisions in 
p+Au, (N co i\) nAu is the average number of collisions in 
n+Au, and dN pp / dq is the yield of the reference nucleon- 
nucleon system. 

The nucleon-nucleon reference comes from the UAl 
measurement m of the p+p inelastic cross section. Note 
that data for p+p is used since data for the preferable 
p+p system is not available at this energy. As described 
in Ref. [B], corrections are applied to the UAl results to 
account for (a) the conversion from rapidity to pseudora¬ 
pidity and (b) the difference between the UAl acceptance 
(\r/\ < 2.5) and the PHOBOS acceptance (0.2 < ij < 1.4). 
An inelastic p+p cross section of 41 mb is used to esti¬ 
mate the yield of p+p collisions given the differential cross 
section measurements from UAl. 

The nuclear modification factor as a function of px in 
the nucleon+Au system, R pn Au , is compared to that of 
d+Au, Rd.Au, for each centrality bin in Fig. [14] Com¬ 
mon systematic errors among the two systems on the 
determination of iV C o]i affect the overall scale of the ra¬ 
tios, as shown by the height of the grey band. Further 
systematic errors in the determination of N co \\ for the 
tagged nucleon+Au system are shown as boxes around 
the Rp n Au points. 

Qualitatively similar results have been found for a nar¬ 
rower window of pseudorapidity by PHENIX [18] , The 
RnAu presented in that work is a simple average of 
p+Au and n+Au, as opposed to the weighted combina¬ 
tion shown in Eq.[4] While the shapes of the modification 
factors are similar in this work and Ref. [18], the latter 
appear to be slightly shifted to larger values, most likely 
due to the use of different reference spectra. 

No significant difference between R pn Au and RdAu 
is observed. This measurement supports the conclu- 
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Figure 14. Comparison of RdAu and R pn Au in each ERing centrality bin. The height of the grey band shows the common scale 
uncertainty due to systematic errors on A T co n. The boxes around the R pn Au points show the supplemental systematic error on 
IVcoii in the nucleon+Au system. 


sions drawn from the nuclear modification factor mea¬ 
surements of g?+Au collisions jS]; namely, that high-px 
hadron production in central Au+Au collisions is signif¬ 
icantly suppressed with respect to the expectation of bi¬ 
nary collision scaling of p+p | 16 j . while the production 
in d+Au collisions is not. It should be noted that no 
claim of binary collision scaling in d+Au or nucleon+Au 
interactions has been made. 

It has been observed that the nuclear modification 
factor in d+Au exhibits a dependence on pseudorapid¬ 
ity mmmm- Thus, the apparent tendency of R pn Au 
and RdAu to take the value of unity at high px is likely a 
consequence of the PHOBOS pseudorapidity acceptance. 
Further, as will be discussed in Sect. VI the hadron pro¬ 
duction of d+Au collisions is known to be enhanced with 
respect to binary collision scaling in a certain range of 
transverse momentum. Any statement that d+Au lacks 
a suppression of high-px hadrons is therefore contingent 
upon the magnitude of this enhancement; see Ref. [49] 
for a discussion. 

Nevertheless, the stark discrepancy observed between 


nucleon+Au and Au+Au collisions at ^/s NN = 200 GeV 
demonstrate that final state effects play a much stronger 
role in the high-px hadron production of central Au+Au 
collisions than do initial state effects. While the pseudo¬ 
rapidity dependence of RdAu may provide evidence of 
some initial modification of the gold nucleus [SB ED, 
it is clear that interactions with some dense, large vol¬ 
ume medium produced only in the nucleus-nucleus sys¬ 
tem forms the dominant source of high-px hadron sup¬ 
pression in Au+Au collisions. The data presented here 
demonstrate that this conclusion is not biased by the use 
of deuteron-nucleus rather than nucleon-nucleus interac¬ 
tions as the control experiment for Au+Au. 


VI. CENTRALITY DEPENDENCE OF THE 
SPECTRAL SHAPE 


Although no clear evidence for enhancements above 
unity are seen in the nuclear modification factor shown 


in Fig. 14 the px dependence may be related to the 
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Figure 15. The ratio of d 2 N c h/dp^dp of d+Au in the 40-70% 
ERing centrality bin to that in p+p, scaled by 2.8 so that the 
ratio is unity at pt = 0.35 GeV/c . The line is the ratio of 
fits to the two spectra using Eq. [l| The arrows mark the pt 
values at which the centrality dependence of the relative yield 
is studied (see Fig. 161. 


so-called Cronin effect. This effect refers to the larger 
ratio of hadron production seen at high px compared 
to lower px in proton-nucleus collisions [7] relative to 
p+p collisions scaled by the effective thickness of the 
nucleus. General aspects of the enhancement of inclu¬ 
sive charged hadron production (that is, unidentified 
hadrons) in p+Au collisions can be described by mod¬ 
els in which partons undergo multiple scattering at the 
initial impact of the p+Au collision |49| . However, the 
observed difference in the strength of enhancement for 
mesons and baryons |52| is not easily explained by ini¬ 
tial state partonic scattering models. While other theo¬ 
ries, such as those based on the recombination model of 
hadronization [53], may be better suited to describe the 
enhancement of individual hadron species, the shape of 
the d+Au px spectrum relative to that of p+p is not 
a thoroughly understood phenomenon. Of particular 
importance is the dependence of the spectral shape on 
the nuclear thickness probed by the projectile (i.e. the 
deuteron in a d+Au collision) [53]. 

The centrality dependence of the nuclear modification 
factor in d+Au and Au+Au collisions at RHIC has been 
studied extensively [30] [55] [57] . A particularly convenient 
method for exploring how the shape of the transverse 
momentum spectra changes relative to p+p has been sug¬ 
gested in Ref. [6]. This method involves studying the cen¬ 
trality dependence of the charged hadron yield in d+Au 
collisions relative to p+p at several values of px- 

The procedure for determining the so-called relative 
yield is as follows. First, the transverse momentum spec¬ 
trum in a particular d+Au centrality bin is compared 


to the spectrum of p+p. To compare only the shape of 
the two spectra, they are then normalized such that the 
spectra match at px = 0.35 GeV/c. While this specific 
value of px is arbitrary, it has been intentionally chosen 
to be in a region where soft processes drive particle pro¬ 
duction. Matching the d+Au spectra to the p+p spectra 
serves to remove any trivial enhancement of hadron pro¬ 
duction in d+Au that is simply due to the larger number 
of nucleon-nucleon collisions occurring in that system. 
However, matching in this way does not assume iV co n 
scaling, nor does it have any effect on the relative shape 
of the spectra. 

Next, the ratio of the normalized d+Au spectra and 
the p+p spectra is determined. The value of this ratio 
at certain transverse momentum values are selected, as 
shown in Fig. [15] Finally, the centrality dependence of 
the normalized ratio, the relative yield, at the chosen px 
values is studied. 


The relative yield of d+Au collisions to p+p is shown 
in Fig. 16 as a function of dN c h/dr), for four different 
values of transverse momentum. It is expected that sys¬ 
tematic effects on the relative yield are highly correlated 
between the spectra measured with different centrality 
bins. Thus, shifts in the relative yield will tend to move 
all points together. See Table [j] for a description of 
the systematic uncertainties on the centrality variables 
measured with ERing. With centrality parametrized by 
the experimentally measured integrated yield, no bias or 
(Glauber) model dependence is introduced by the choice 
of centrality variable. 


From Fig. m it is clear that the difference between the 
d+Au and p+p spectra depends on both centrality and 
px- If the shape of the two spectra were identical, the 
relative yield would be constant at unity for all values 
of px and centrality. Instead, the d+Au spectra show 
an enhancement over p+p that increases with centrality. 
The strength of this enhancement is observed to increase 
at higher px- It would be interesting to study the rel¬ 
ative yield of much higher px hadrons, on the order of 
10 to 100 GeV/c, in order to test whether the shape of 
the p+p spectra is recovered in hard scattering processes. 
However, such particles are produced very rarely and too 
few are present in the PHOBOS data set to allow such a 
study. 

Nevertheless, the data show a smooth extrapolation 
of the relative yield of d+Au collisions to that of p+p 
as the d+Au collisions become more peripheral. Thus, 
distortions of the d+Au spectra caused by nuclear ef¬ 
fects diminish in a smooth way as the amount of nuclear 
material probed by the deuteron is reduced. The inte¬ 
grated charged particle yield near (p) ss 0.8 has been cho¬ 
sen as the centrality measure, since it provides a model- 
independent variable with which to study the centrality 
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Figure 16. Open symbols show the average hadron yield of d+Au collisions relative to p+p as a function of d+Au dN^/dp 
near mid-rapidity, scaled by a factor, fc, such that the ratio is unity at pr = 0.35 GeV/c in order to focus on the evolution of 
the shape of the yield. Statistical errors are represented by bars on the points. The systematic error for the ratio in ERing 
centrality bins is shown by the grey band. See text for a discussion of the systematic errors. Closed triangles at a relative yield 
of 1 (representing p+p divided by itself) are plotted at the d/V c h /dp for p+p. The dependence of the relative yield on both 
centrality and pt is observed to extrapolate smoothly back to p+p. 


dependence of hadron production in nucleon-nucleus and 
nucleus-nucleus systems. 


VII. CHARGE TRANSPORT 

The availability of both p+Au and n+Au collision data 
presents a unique opportunity to study baryon transport 
in nucleon-nucleus collisions. Since a p+Au collision con¬ 
tains one more charged hadron than an n+Au collision, 
a search for this extra charge near the mid-rapidity re¬ 
gion is possible. Previous measurements [35] of p+Au 
collisions at ^ NN = 19.4 GeV found that the number 
of net protons (p - p) per unit of rapidity is less than 
one in the mid-rapidity region. In addition, studies have 
shown a decrease in the mid-rapidity net proton yield 


with increasing center of mass energy; see Ref. [55] for 
a discussion. Further, it has been inferred that hadrons 
traversing nuclear material do not lose more than about 
two units of rapidity [60]. Thus, it is expected that any 
charge asymmetry between hadrons measured at mid¬ 
rapidity in p+Au and n+Au collisions would be small. 

Nevertheless, a comparison of charged hadron pro¬ 
duction in p+Au and n+Au allows the transport of 
charge explicitly from the projectile proton to be stud¬ 
ied. Assuming that baryons from the gold nucleus un¬ 
dergo transport to mid-rapidity via the same process in 
both p+Au and n+Au collisions, any charge transport 
to mid-rapidity of protons in the gold nucleus would not 
lead to an asymmetry. 

Simple charge conservation would imply that the total 
number of positive particles emerging from a p+Au colli- 
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sion should be greater (by one) than the number emerg¬ 
ing from a n+Au collision. Whether or not this charge 
asymmetry is present near mid-rapidity has been studied 
using the observable , defined as 

AP n = dNlt'/dn-dN^/dn 
h± dNli u / dq + dN^ u / dq 

where A v } ^_ denotes the asymmetry between p+Au and 
n+Au in the yield of positively charged hadrons at 
( 77 ) = 0.8 and A v } ^_ denotes the asymmetry of the yield of 
negatively charged hadrons between the two systems. 
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Figure 17. The asymmetry of positive hadrons between 
p+Au and n+Au collisions at (rj) = 0.8 as a function of cen¬ 
trality. The grey band shows the systematic uncertainty in 
the overall scale of the ratio. 
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Figure 18. The asymmetry of negative hadrons between 
p+Au and n+Au collisions at (rj) = 0.8 as a function of cen¬ 
trality. The grey band shows the systematic uncertainty in 
the overall scale of the ratio. 


The charge asymmetry defined by Eq. [5| is presented 
in Fig. [XT] for positive hadrons and in Fig. [18] for negative 
hadrons. The grey band in each figure represents the sys¬ 
tematic uncertainty in the asymmetry ratio, propagated 
from the nucleon tagging component of the systematic 
uncertainty on the momentum spectra (see Sect. IVI. 
Only uncertainties specific to reconstructing the nucleon- 
nucleus px spectra contribute to this systematic error, as 
all other effects divide out in the ratio. No evidence for 
asymmetry between p+Au and n+Au collisions is ob¬ 
served at (p) = 0 . 8 , which is slightly forward on the 
deuteron-going side. 


VIII. SUMMARY 


The addition of two forward proton calorimeters to 
the PHOBOS detector allows the extraction of p+Au 
and n+Au collisions from the d+Au data set. Centrality 
parameters have been determined for each of the colli¬ 
sion systems using observables based on the multiplic¬ 
ity at high rapidity and on the number of spectators. 
The number of particles produced near mid-rapidity is 
found to scale with N ps , rt across all collision systems. 
The charged hadron spectra have been measured for 
p+Au, n+Au, and d+Au collisions and used to construct 
an ideal nucleon-nucleus reference for Au+Au collisions. 
The nuclear modification factor of this ideal reference is 
found to agree with that of d+Au. The shape of the nu¬ 
clear modification factor has been studied in detail and 
is found to depend on both centrality and transverse mo¬ 
mentum. A larger ratio of the d+Au over p+p spectra 
is found at larger values of px and this enhancement is 
found to extrapolate smoothly as a function of multiplic¬ 
ity at mid-rapidity from p+p to central d+Au collisions. 
Finally, a comparison of the yield of positively and neg¬ 
atively charged hadrons in p+Au and n+Au has been 
conducted in a direct search for evidence of charge trans¬ 
port to mid-rapidity. No significant asymmetry between 
the charged hadron yields in p+Au and n+Au is observed 
at ( 77 ) = 0 . 8 . 
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